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ABSTRACT 

Context. Observations of the cosmic microwave background, light element abundances, large-scale distribution of galaxies, and dis- 
tant supemovae are the primary tools for determining the cosmological parameters that define the global structure of the Universe. 
Aims. Here we illustrate how the combination of observations related to strong gravitational lensing and stellar dynamics in elliptical 
galaxies offers a simple and promising way to measure the cosmological matter and dark-energy density parameters. 
Methods. A gravitational lensing estimate of the mass enclosed inside the Einstein circle can be obtained by measuring the Einstein 
angle, once the critical density of the system is known. A model-dependent dynamical estimate of this mass can also be obtained by 
measuring the central velocity dispersion of the stellar component. By assuming the well-tested homologous 1/r- (isothermal) profile 
for the total (luminous+dark) density distribution in elliptical galaxies acting as lenses, these two mass measurements can be properly 
compared. Thus, a relation between the Einstein angle and the central stellar velocity dispersion is derived, and the cosmological 
matter and the dark-energy density parameters can be estimated from this. 

Results. We determined the accuracy of the cosmological parameter estimates by means of simulations that include realistic mea- 
surement uncertainties on the relevant quantities. Interestingly, the expected constraints on the cosmological parameter plane are 
complementary to those coming from other observational techniques. Then, we applied the method to the recent data sets of the Sloan 
Lens ACS (SLAGS) and the Lenses Structure and Dynamics (LSD) Surveys, and showed that the concordance value between 0.7 and 
0.8 for the dark-energy density parameter is included in our 99% confidence regions. 

Conclusions. The small number of lenses available to date prevents us from precisely determining the cosmological parameters, but 
it still proves the feasibility of the method. When applied to samples made of hundreds of lenses that are expected to become avail- 
able from forthcoming deep and wide surveys, this technique will be an important alternative tool for measuring the geometry of the 
Universe. 

Key words, cosmology: theory - cosmology: observations - galaxies: distances and redshifts - galaxies: elliptical and lenticular, cD 
- gravitational lensing - galaxies: kinematics and dynamics 



1. Introduction Ho = (72 + 8)kms"' Mpc"' (Freedman et al.lSOQll). The esti- 

\^ ' mates of different cosmological parameters from a single obser- 

C^ . The Universe appears to be dominated by dark-energy and dark ^^^^^^^^j ^^^^^^ ^^.^ ^^^^^ con-elated (hence "degenerate") and 

■ ■ matter. Although the physical nature of these dark components ^^^^-^ significant uncertainties. For instance, from the WMAP 

IS still unknown, the standard cosmological ACDM model with ^j^^^^ ^^^^ ^j^^^^ ^^^1^^^^^ ^ -^^ ^^ ^j^^ ^^^^^^^ ^f jj^^ 

only a few parameters fits most of the current data well: preci- Universe, the best-fit model is characterized by Q„, = 0.42, 

sion measurenients of the anisotropiesm the cosmic microww ^^ ^ 0^3^ Hq = 55 km s' Mpc' (Spergel et al. EQQ3, val- 

background (CMB; Bennett et al. ^0031 Spergel et al. |2003l ^^^ ^j^^^ ^^.^ -^^ ^^^^^^^ ^^^^ ^^e concordance values reported 

^, the observed abundances of hght elements (Buries et al. ^^^^^ ^^-^ ^^^^ -^^ measurements of the cosmo- 

EQpTi Cyburt et al. 2003 ) the ai-ge-scale distribution of galaxies j ^^^j g^rs can only be obtained by using complemen- 

(LSS; Tegmark et al. 2004 , Cole et al. 2005 ), and the luminosity- techniques. In fact, considerable efforts are still being made 

distance relationship for distant typela supernovae (SNIa; Riess -^ ^^^^^ ^^ ^^^^^^ ^^^^^^^^ measurements of these parameters (in 

et al [SmEOQll Perlmutter et al ^1999^. In this standard model, ^1,^1,^, see the scientific goals of the forthcoming PLANCK 

the Universe is homogeneous and isotropic on its lai'gest scales, . cam n missions) 
and its geometry appears to be flat (Q = Q„, -1- Qa ~ 1); 

the total mass-energy density is mainly in the form of dai'k- The deflection of light due to gravitational lensing is sensi- 

energy (Qa * 0.7) and matter (Q„, ^ 0.3), ordinary and dai'k. tive to the total matter density of the structures in the Universe, 

These values of the cosmological parameters imply a fairly re- independendy of the nature or dynamical state of the deflecting 

cent transition from a decelerating to an accelerating univer- mass. Therefore, strong and weak gravitational lensing provide 

sal expansion. The Hubble Space Telescope (HST) Key Project valuable tools for measuring the distribution of mass. In partic- 

has measured the current expansion rate, the Hubble parameter ular, cosmic shear estimates of the amplitude of the weak lens- 
ing distorsions of distant sources over a wide range of angular 



Send offprint requests to: C. Grillo scales (Bartelmann & Schneider 1 19991 ) are a very encouraging 



C. Grillo et al.: Cosmological parameters from strong gravitational lensing and stellar dynamics in elliptical galaxies 



way to study the large-scale structure of the Universe, and there- 
fore to probe the parameters that define the relevant cosmolog- 
ical model (Refregier 2003). Further constraints on the geom- 
etry of the Universe can also be provided by the abundance of 
lenses or arcs observed in lens surveys (see Bartelmann & Weiss 
119941 for a numerical approach; Mitchell et al. 120051 for obser- 
vational results from the CLASS Survey described by Myers et 
al. [2003 and Browne et al. ^003 ) and by strong (Link & Pierce 
[19981 Soucail et al. [2004b or weak (Lombardi & Bertin [19991 
Jain & Taylor I2003I I lensing mass reconstructions in clusters 
of galaxies. Isolated strong lens galaxies offer another possi- 
bility to measure the cosmological parameters (Kochanek |1992[ 
11996; Myungshin et al. 1997) . including the Hubble parameter 
(Refs dal [T9M1 Koopmans et al. I2003bt Mortsell & Sunesson 
[20061 ). 

The various techniques that have been proposed are limited 
by several assumptions. For example, a measurement of the mat- 
ter density parameter from cosmic shear is degenerate with that 
of the normalization of the amplitude of the power spectrum of 
matter perturbations (erg); moreover, an extremely large number 
of high-quality galaxy images and some modeling on the growth 
of the structure in the Universe are required. For techniques 
based on gravitational lens statistics, the luminosity function, 
the relation between luminosity and velocity dispersion, and the 
density profile of the lens galaxies play an important role. The 
use of arcs statistics awaits more realistic simulations of clusters 
and observations, since different studies have led to contrasting 
results (e.g., Bartelmann et al. |1998l Dalai et al., 2004j . The esti- 
mates of the cosmological parameters from cluster mass recon- 
structions have also two distinct limitations: the presence of pos- 
sible substructure in the region of multiple image formation (for 
strong lensing) and the need for a high density of background 
galaxies (for weak lensing). 

In this paper we propose a new technique that, starting from 
strong gravitational lensing and stellar dynamics observations in 
elliptical galaxies, is able to probe the geometry of the Universe 
in a different and effective way. The paper is organized as fol- 
lows. In Sect. 2, we describe our method to estimate the matter 
and dark-energy density parameters. Then, in Sect. 3, we deter- 
mine through simulations the precision attainable in these mea- 
surements of the cosmological parameters. Good estimators for 
the quantities relevant to the problem are identified in Sect. 4. 
In Sect. 5, we apply our technique to the data collected and just 
published of two surveys of lens galaxies for which stellar dy- 
namical measurements are available. Finally, in Sect. 6, we sum- 
marize and discuss the results obtained. 

2. The method 

It is known that in an axisymmetric lens multiple images can 
only form in the vicinity of the so-called Einstein ring, at an an- 
gle 0E from the center of the lens (see Schneider et al. 119921 ). 
From the theory of gravitational lensing, the "mass" ./#gri en- 
closed within the disk defined by the Einstein ring is directly re- 
lated to the geometry of the configuration, through the definition 
of the critical density (^gri - 2ci-7r0g). This "mass" ^g,-i is con- 
nected to the intrinsic mass Mgri of the lens by the distance to the 
lens (by converting the Einstein angle into an Einstein radius). 

A dynamical estimate ^dyn of the mass can also be given 
by measuring the quantity (TqOe, where (Tq is the central velocity 
dispersion of the stellar component, usually referred to the disk 
of radius /?e/8 (Re being the standard optical effective radius); 
the dynamical mass is then obtained by multiplication by a suit- 
able factor (^dyn - o cTqOe) that is model-dependent (the lens 



usually includes a significant dark matter component). Again, in 
order to relate ^dyn to the intrinsic mass Mdyn we should convert 
9e into a radius. 

With no need to refer to intrinsic masses (and thus with no 
need to know the exact distance to the lens, which would bring 
in a knowledge of the Hubble constant //q), we thus see that, if 
we identify ^g^i = ^dyn, the combination of a measurement 
of 0E and of ctq should be uniquely related, in the standard cos- 
mological model, to a function of the redshifts zi and Zs (of the 
lens and of the source, respectively) and of the cosmological pa- 
rameters Q,„ and Q.^- Note that, in this context, zi and Zs can be 
considered to be measured with negligible errors. 

Several studies, based on various dynamical tracers (stars, 
globular clusters, planetary nebulae. X-ray halos, HI disks and 
rings), have established that bright elliptical galaxies of the lo- 
cal universe, as a rule, exhibit approximately flat circular ve- 
locity curves (e.g., Gerhard et al. 2001, Peng et al. 2004), thus 
suggesting that the structure of these systems should be consid- 
ered as approximately homologous, with the total density dis- 
tribution (luminous-i-dark) close to that of a singular isothermal 
sphere (SIS; p oc l/p-). These detailed dynamical studies refer 
to nearby galaxies and thus the result obtained does not depend 
on the values of the cosmological parameters. Of course, this 
is a zeroth-order description, and different galaxies may exhibit 
different deviations from this "universal" total density profile. 

A number of investigations of galaxies at cosmologically sig- 
nificant distances address observed properties that result from 
the combined effects of evolution and of the geometry of the 
universe. The interpretation of these data can thus be obtained 
in different ways. For example, in the study of the Fundamental 
Plane out to z «; 1 (e.g., see Treu et al. 2002, and other parallel in- 
vestigations, many of which quoted there), one may assume ap- 
proximate structural homology and a given cosmological model 
and use the data on the observed change in the Fundamental 
Plane to derive information on the evolution properties of the 
observed stellar populations. 

Here we recall that, under the assumption of the concor- 
dance cosmological model (Q„, - 0.3 , Qa = 0.7 , Hq = 
70 km s"' Mpc"^), analyses of strong gravitational lensing alone 
(e.g., Rusin e al. 120031 ) or combined measurements of stellar dy- 
namics and gravitational lensing (e.g., Treu & Koopmans |2004t 
Koopmans et al. 120061 ) have confirmed the persistence of struc- 
tural homology, i.e., that little evolution appears to take place 
in the observed total density profile of bright ellipticals, in the 
sense that they are found to be characterized by approximately 
SIS density profiles also at cosmologically significant distances. 
In Sect. 5 we will show that this conclusion is robust with respect 
to the choice of the adopted cosmological parameters. 

These findings have encouraged us to explore the conse- 
quences of considering the combined measurements of stellar 
dynamics and gravitational lensing on distant ellipticals, start- 
ing from the simplifying assumption that homology is indeed 
strictly followed by these systems; in particular, we wish to ex- 
plore whether this assumption, applied to the interpretation of 
the data, may lead to interesting constraints on the values of the 
parameters that define the geometry of the universe. In other 
words, we wish to check the consequences on the cosmologi- 
cal parameters of assuming from the very beginning that (1) the 
espression for ^gri is basically that associated with an SIS, and 
(2) no significant variation on the "virial coefficient" a is present 
from galaxy to galaxy. 
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In practice, we proceed as follows. We note that for an SIS 



Distance ratio (£!,„ = 0.3, £J^ = 0.7; z, var.) 



0E = 47r 



o"sis\ Djs_ 

Dos 



(1) 



where crgis is the lens "velocity dispersion", c is the light speed, 
Dis and Dgs are the lens-source and the observer-source angu- 
lar diameter distances, respectively. We then bypass the issues 
raised by dynamical modeling by recalling that ctq turns out to be 
a good estimate of osis . This latter point, exploited by Kochanek 
(11993 1994), was confirmed by Treu et al. (2006), and is now 
checked by us independently, by a test described separately in 
Sect. 4.2 on a sample of eight well-studied nearby ellipticals. 
Therefore we consider the quantity 



47r 



1^ 



Dh_ 
Do. 



= r(z/,z^;Q,„,QA), 



(2) 



as the observable that will be used to produce, by studying a 
statistically significant sample of lenses, a measurement of fi,„ 
and i^A. 

Given the weak dependence of the relevant coeflicients, such 
as a, on the detailed relative distributions of dark and lumi- 
nous matter in galaxies, the method is expected to be robust. 
In particular, the entire argument could be easily generalized to 
non-axisymmetric lenses by referring to the properties of the so- 
called singular isothermal ellipsoid (see Kormann et al. |19941 l. 

Figure[T]illustrates the distance ratio r of Eq. (|2|i (which does 
not have a simple analytic form) versus the source redshift, as a 
function of the lens redshift and of the matter and dark-energy 
density parameters. In general, different cosmological models 
give values of r which differ more clearly at higher values of 
the source redshift (see the last three panels); in addition, the 
quantity r is more sensitive to small variations of Q.j^ than of fi,„ 
(compare the second and third panels). As a consequence, we 
naturally expect this method to be optimally efficient in measur- 
ing Qa, provided we have at our disposal a sample of lenses at 
sufficiently high redshifts. 
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Fig. 1. Dependence of the angular diameter distance ratio r = 
D Is I Dos on the lens redshift {top left), matter density parame- 
ter (top right), dark-energy density parameter {bottom left), and 
matter and dark-energy density parameters in a flat cosmologi- 
cal model {bottom right). In the first panel, (Q,„, Qa) - (0.3, 0.7) 
are fixed and zi is varied from 0.1 to 1, with a regular step of 
0.1. The arrow shows the increasing direction of z;. In the sec- 
ond panel, Qa = 0.7 and zi - 0.5 are fixed and Q,„ is varied from 
to 1 , with a regular step of 0. 1 . The arrow shows the increasing 
direction of Q„,. In the third panel, Q.„i - 0.3 and zi - 0.5 are 
fixed and Qa is varied from to 1 , with a regular step of 0. 1 . The 
arrow shows the increasing direction of Qa- In the fourth panel, 
Q„, + Qa = 1 and zi = 0.5 are fixed and Qa is varied from 
to 1, with a regular step of 0.1. The arrow shows the increasing 
direction of Qa- 



3. Simulated measurements of the cosmological 
parameters 

In order to explore the precision with which the above technique 
can probe the cosmological parameters, we have performed sev- 
eral simulations. We modelled each lens as a singular isothermal 
sphere with an external shear component. Because of the known 
degeneracy between external shear and ellipticity (Witt & Mao 
119971 1. we remark that our modeling choice is also good at de- 
scribing simulations with singular isothermal ellipsoid models. 
As suggested by real lensing systems, we considered the lens 
redshift uniformly distributed between and 1 and the source 
redshift between 1 and 3.5; the lens velocity dispersion and the 
external shear values were drawn from uniform distributions 
ranging from 100 to 350 kms ' for the first, from to 0.2 in 
magnitude and from 0° to 180° in orientation for the second (the 
role of the external shear component will become relevant in 
Sect. 4.1). We calculated the Einstein angles from Eq. ([T) in a 
(QotjQa) - (0.3,0.7) cosmology and completed each lensing 
system simulating randomly, inside a square of side 6e and cen- 
tered on the lens, the position of a source. Only the lenses with 
two images of this source and with 6e greater than 0.5" were 
accepted for the next analyses. The reason is that we wanted 
to investigate multiple image systems similar to those observed. 



where the images are far enough from the center of the galaxy 
acting as a lens to be resolved with the present technology. 

We started by examining the dependence of the error es- 
timates on Qm and Qa on the simulated uncertainties on the 
Einstein angle and on the central velocity dispersion. In order 
to do this, we minimized a chi-square {x^) like estimator with 
respect to the two cosmological parameters. This function is de- 
fined by comparing the "observational" distance ratio calculated 
from 6e and o-q, that is the quantity on the left in Eq. dU, and the 
"theoretical" distance ratio obtained from the lens and source 
redshifts, that is the quantity on the right in Eq. (|2]i: 



X 



-(Q,„,Qa) = ^ 



fjp^ -r(z/,,,Zj,;Q„„QA) 



"(fi)(t)<*,)^-fe)(K)" 



(3) 



We included realistic measurement errors on the Einstein angle 
{SOe) and on the central velocity dispersion {dcro), considering 
normal distributions with standard deviations equal to 4%, 5%, 
6%, and 7% of the true values. At first, we performed 2000 min- 
imizations for samples of 100 and 200 lenses, assigning a nom- 
inal 0% uncertainty to one of the two quantities and varying the 
error on the other in the range reported above. Then, we repeated 
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Parameter distribution (4% err. on 9^) Parameter distribution (4% err. on (J^) Parameter distribution (4% err. on 9^) Parameter distribution (4% err. on <Jq) 
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Fig. 2. Estimates of the cosmological parameters. Simulation of 2000 measurements composed of 100 {on the left) and 200 {on 
the right) lenses each, with different uncertainties (increasing from the top to the bottom) on the Einstein angle {first and third 
column) and on the central velocity dispersion {second and fourth column). A nominal 0% uncertainty is assigned to the quantity 
not mentioned in the panels. Thick bars on the co-ordinate axes and contour levels on the planes represent, respectively, the 95% 
confidence intervals and the 68% and 95% confidence regions for the cosmological parameters. A cross shows the position of the 
true parameters: (Q„„ Qa) = (0.3, 0.7). 



the analysis just described with the additional hypothesis of a flat which these quantities can be measured at present by means of 
cosmological model (Q,„ + Qa = 1)- Finally, we simulated es- the best ground and space-based telescopes, 
timates of the cosmological parameters, in general and flat cos- 
mology models, considering 5% eiTors on both the Einstein an- The results for the joint and single probability density func- 
gle and the central velocity dispersion. This is the precision with tions (/) of the cosmological parameters are summarized by 

Figs.|2l[3] and Hi and by Tables [T]|2] and [3] As a first good in- 
dication, the minimum jt'^ was found to be asymptotically unbi- 
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Fig. 3. Estimates of the cosmological parameters in a flat model (Q.^ + Qa = 1). Simulation of 2000 measurements composed of 
100 (on the left) and 200 {on the right) lenses each, with different uncertainties (increasing from the top to the bottom) on the 
Einstein angle (first and third column) and on the central velocity dispersion (second and fourth column). A nominal 0% uncertainty 
is assigned to the quantity not mentioned in the panels. Thick bars on the abscissa axes represent the 95% confidence intervals for 
the cosmological parameters. Crosses show the position of the true parameters: (Q.m, ^^a) - (0.3, 0.7). 



ased, i.e. centered on the original values (Q.m, ^\) - (0.3, 0.7) in 
the limit of small uncertainties on 9^ and ctq. Then, we remark 
that the true values of the cosmological parameters are always 
inside the regions and intervals at the 95% confidence level, but 
their boundaries are not symmetrical. In fact, in all the three pre- 
viously mentioned tables these confidence intervals for Q„, are 
more extended on the high side; in contrast, the intervals for Qa 



are generally larger on the low side. Moreover, since the dis- 
tance ratio r depends linearly on the Einstein angle and quadrati- 
cally on the central velocity dispersion, larger confidence regions 
and intervals are obtained when a fixed uncertainty is considered 
on the latter quantity. The assumption about the flatness of the 
Universe makes our technique more powerful, allowing a better 
estimate of the dark-energy density parameter, as can be seen 
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Table 1. Intervals at 95% confidence level for the matter and the dark-energy density parameters. 
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5% 


6% 


7% 


Err. (To 


4% 


5% 


6% 


7% 


^A 


[0.00,0.81] 
[0.58,0.78] 


[0.00,0.98] 
[0.54, 0.80] 


[0.00,1.19] 
[0.50, 0.80] 


[0.00, 1.41] 
[0.47,0.80] 


Ha 


[0.03,0.63] 
[0.60, 0.75] 


[0.00,0.73] 
[0.57,0.75] 


[0.00,0.85] 
[0.54, 0.75] 


[0.00,0.97] 
[0.50,0.74] 



Notes - These intervals are obtained by excluding from the 2000 ;t'^ minimizations the 50 smallest and the 50 largest values for fi,,, and Q.a- The 
simulated measurement uncertainties of the Einstein angle (on the top) and of the central velocity dispersion (on the bottom) range from 4% 
to 7% (a nominal 0% uncertainty is assigned to one of the two quantities and the error on the other is varied) for samples of 100 (on the left) 
and 200 (on the right) lenses. 

Table 2. Intervals at 95% confidence level for the matter and the dark-energy density parameters in a flat cosmology (0,„ + Qa = 1). 







N=100 










N=200 






Err. ^E 


4% 


5% 


6% 


7% 


Err. ^E 


4% 


5% 


6% 


7% 


^A 


[0.28,0.34] 
[0.66, 0.72] 


[0.28,0.36] 
[0.64, 0.72] 


[0.28,0.38] 
[0.62, 0.72] 


[0.28,0.40] 
[0.60, 0.72] 


Ha 


[0.29,0.33] 
[0.67,0.71] 


[0.29,0.34] 
[0.66,0.71] 


[0.29,0.36] 
[0.64,0.71] 


[0.30,0.38] 
[0.62, 0.70] 


Err. ao 


4% 


5% 


6% 


7% 


Err. (To 


4% 


5% 


6% 


7% 


Ha 


[0.26,0.38] 
[0.62, 0.74] 


[0.26,0.41] 
[0.59, 0.74] 


[0.26,0.44] 
[0.56, 0.74] 


[0.26,0.49] 
[0.51,0.74] 


Ha 


[0.28,0.36] 
[0.64, 0.72] 


[0.28,0.39] 
[0.61,0.72] 


[0.28,0.41] 
[0.59,0.72] 


[0.29,0.45] 
[0.55,0.71] 



Notes - These intervals are obtained by excluding from the 2000 ;t' minimizations the 50 smallest and the 50 largest values for Q.„, and Q./^. The 
simulated measurement uncertainties of the Einstein angle (on the top) and of the central velocity dispersion (on the bottom) range from 4% 
to 7% (a nominal 0% uncertainty is assigned to one of the two quantities and the error on the other is varied) for samples of 100 (on the left) 
and 200 (on the right) lenses. 
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Fig. 4. Estimates of the cosmological parameters assuming 5% 
uncertainty on both the Einstein angle and the central veloc- 
ity dispersion. Simulation of 2000 measurements composed of 
100 (on the top) and 200 (on the bottom) lenses each, in gen- 
eral (on the left) and in flat (on the right) cosmological mod- 
els. Thick bars on the co-ordinate axes and contour levels on 
the planes represent, respectively, the 95% confidence intervals 
and the 68% and 95% confidence regions for the cosmological 
parameters. Crosses show the position of the true parameters: 
(Q,„,Qa) = (0.3,0.7). 



Table 3. Intervals at 95% confidence level for the matter and 
the dark-energy density parameters in general (second and third 
columns) and flat (fourth and fifth columns) cosmology. 



ii,,j 



n. 



n,„ 



Q.A- 



N= 
N= 



100 
200 



[0.00,1.12] 
[0.00,0.80] 



[0.52,0.79] 
[0.54,0.74] 



[0.27, 0.45] 
[0.29,0.42] 



[0.55,0.73] 
[0.58,0.71] 



" Flat cosmology: Q.„ + Ha = 1- 
Notes - These intervals are obtained by excluding from the 2000 x~ 
minimizations the 50 smallest and the 50 largest values for Q.„, and 
Da- The simulated measurement uncertainties are 5% of the true 
values for both the Einstein angle and the central velocity disper- 
sion. Samples of 100 and 200 lenses are considered. 



by comparing Tables [T] and |2] Finally, from the results of Fig. 
|4]and Table [3] we argue that the cuiTent measurement precision 
akeady allows a good estimate of Qa from a sufficiently large 
sample of lenses. 

From what discussed so far, we conclude that this method 
is especially well-suited to measure the dark-energy density pa- 
rameter. This is not surprising, given the comments made in the 
previous section about the different dependence of the distance 
ratio on the two cosmological parameters. In view of the fact 
that, so far, only supernova observations have been able to reveal 
directly a significant dark-energy component, our technique of- 
fers a promising new test for the concordance model. It is also 
particularly interesting to notice that the 95% confidence regions 
on the parameter planes in Fig. |2] and Fig. |4] are oriented in such 
a way to be complementary to the results of the other cosmolog- 
ical probes currently considered. 
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Fig. 5. Consistency and bias tests for the 0e estimators. Top: 
Angular coefficients (mi and 1112) of the linear fits of the true 
Einstein angle versus the estimated value from 0i (on the left) 
and 02 {on the right), for an increasing number of lenses. Bottom: 
Angular coefficients of the linear fits of the true Einstein angle 
versus the estimated value from 0i {on the left) and 02 {on the 
right) for 10000 samples with 100 lenses each. Solid lines show 
the mean values. 

Table 4. Statistics of the two estimators of the Einstein angle. 



Bias 



MSE 



6*1 

02 



-0.020 
-0.008 



0.0008 
0.0005 



4. Diagnostics of the relevant quantities 

4. 1 . The Einstein angle 

Isolated lens galaxies which show low ellipticity values in the 
luminous distribution are expected to be well described by ax- 
isymmetric models; on the other hand, non-axisymmetric lens 
models are required to represent galaxies in groups or clusters 
and galaxies displaying high ellipticity in the luminosity profile. 
In the former case, the lens properties are embodied by the so- 
called Einstein angle; in the latter case, an equivalent Einstein 
angle can still be considered. However, in both cases the un- 
certainty with which the Einstein angle can be reconstructed by 
lensing modeling is on the order of 5%, provided high-quality 
imaging of the systems is available. 

In future cosmological studies we will have to handle data 
sets from wide and deep sky surveys made of hundreds of lenses. 
In order to evaluate in a simple way 0e for lensing systems with 
two images of the same source, we may define two estimators. 
They could be used to perform a preliminary fast analysis of the 
data, before detailed models with better precision are built. (This 
problem has not been taken into account in the measurements of 
the cosmological parameters presented in this paper, because in 
the following we will consider only data for which refined lens 
models are available). The first estimator {0i) is given by the 
semi-distance between the two images, and the second one (62) 
by the semi-sum of the distances of the two images from the cen- 



ter of the lens. The triangular inequality ensures that 02 is always 
greater than 0i . We remark that the first quantity is easy to mea- 
sure, even when the lens (mass) center is not known a priori; any- 
way, in most cases the lens center should be identified with the 
galaxy luminosity centroid. The two estimators coincide when a 
single axisymmetric model is considered, but they differ when 
some ellipticity for the lens or an external shear component is 
present. In Fig.|5]and Table|4]we illustrate the statistical proper- 
ties of 01 and 02, obtained from simulations in which the physical 
variables were selected as described in the previous section. The 
consistency test and the values of bias and mean squared error 
(MSE) (for definitions see Cowan 1 19981 1 favour 02 as the better 
estimator for 0e- 



4.2. Tine velocity dispersion 

The velocity dispersion of stars in galaxies is a well-defined dy- 
namical quantity (see Bertin '2000'). The expression "velocity 
dispersion" is also used in lensing studies to refer to a param- 
eter of the 1 /r^ density distribution characteristic of the isother- 
mal sphere. In fact, the measured total mass of a lens within the 
Einstein angle can be easily converted into an effective velocity 
dispersion, given the relation between 0e and crsis of Eq. ([T]i- As 
already mentioned in Sect. 2, a one-component isothermal model 
has proved to be an adequate description for the total density 
distribution in elliptical galaxies, as far as lensing is concerned. 
On the other hand, dynamical modeling requires two-component 
(luminous+dark) models in order to find valid agreement with 
the observations; in principle, the velocity dispersion of the stel- 
lar component might be different from that of the dark compo- 
nent. Here we wish to compare the central value of the stellar ve- 
locity dispersion to the velocity dispersion of a one-component 
isothermal model supposed to represent the total (stellar and 
dark) matter distribution. We will find that a good estimator of 
the one-component velocity dispersion of an isothermal model 
(csis) is indeed the stellar central velocity dispersion {o-q). This 
step is essential in order to apply the method proposed in this 
paper. 

We have considered eight bright, nearly round, early-type 
galaxies, which were modeled as described in Bertin et al. ( 119921 ) 
and Saglia et al. (1992'). We assumed that the total mass of 
each galaxy could be described in terms of a singular isother- 
mal sphere and calculated crsis starting from the best-fit two- 
component models of Saglia et al. ( 1992). Then, using the kine- 
matical data reported in Davies & Birkinshaw ( 11988 ) and Franx 
et al. ( 119891 ), we estimated o-q for the same galaxy sample. The 
results are summarized in Table |5j the ratio of the two velocities 
{q - CTo/cTsis) is displayed in Fig.|6] The uncertainty on ctq is just 
the rms scatter of the velocity dispersions measured at different 
slit position angles, so it is probably an underestimate of the real 
value. In addition to this, a minimum 5% error on the crsis was 
assumed in order to obtain the error bars of Fig.|6] We notice that 
q has an average value very close to unity (1.003 + 0.017) and 
an intrinsic rms scatter remarkably small (0.047). Therefore, we 
infer that ctq, which is related to the stellar component alone, is 
a good diagnostic of crsis, a tracer of the total (stellar and dark) 
matter distribution. Moreover, the ratio q does not show any sig- 
nificant dependence either on the crsis value, as can be seen in 
FiglH or on the dark matter fraction /dm inside 7?e, as reported 
in Table |5] This fact suggests the existence of an efficient mech- 
anism of coupling between stellar and dark mass ("conspiracy"). 
In addition, we also remark that q is not influenced by a partic- 
ular choice of the cosmological parameter values, because the 
galaxy set is placed in the nearby Universe. 
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Table 5. Velocity dispersions and dark matter fraction for a sam- 
ple of eight bright, nearly round, early-type galaxies. 



Table 6. The lens samples of the SLACS (1-15) and LSD (16- 
20) Surveys. 



Object 


NGC 1404 


NGC 1549 


NGC 3379 


NGC 4278 


o-sis 
/dm 


245 ±7 
259 
0.58 


198 ±6 
196 
0.59 


219 ± 13 
212 
0.40 


228 ±9 
232 
0.18 


Object 


NGC 4374 


NGC 4472 


NGC 4486 


NGC 4636 


o"sis 
/dm 


291 ±6 
280 
0.59 


292 ±7 
302 
0.33 


311 ±3 
286 
0.34 


206 ±8 
213 
0.27 



Notes - Stellar central velocity dispersion ((Tq), one-component veloc- 
ity dispersion of the singular isothermal sphere (iTsis) that best fits 
the photometry and the kinematics, and dark matter mass fraction 
inside R^ (/dm). 

References - Davies & Birkinshaw ( 11988b : Franx et al. ( il989D : Saglia 
et al. ( fT992l >. 
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Fig. 6. Ratio (q) of the stellar central velocity dispersion (o-q) to 
the one-component velocity dispersion of the singular isother- 
mal sphere (crsig) that best fits the photometry and the kinemat- 
ics for a sample of eight bright, nearly round, early-type galax- 
ies. 



A larger sample (the SLACS sample, described in the fol- 
lowing section) of more distant early-type galaxies, but with es- 
sentially the same range of velocity dispersion, was studied sim- 
ilarly by Treu et al. (.2006) . There they estimated the isothermal 
velocity dispersion inside the Einstein radius Re from the best 
lensing model and their results about the value and intrinsic scat- 
ter of q are consistent with what we found. In conclusion, since 
we demonstrated that cro is a good estimator of crgis, we argue 
that the former quantity can be substituted in Eq. ([TJ to get the 
cosmology-dependent relation given by Eq. (|2]l. 

5. The SLACS and LSD samples 

As an application of our technique, we considered the SLACS 
(Sloan Lens ACS) Survey data set presented in Bolton et al. 
dlOOel l. Treu et al. (|2006] |. and Koopmans et al. (|2006] | (see 
Table|6]l. The sample is composed of fifteen massive, early-type. 



# 



Ts e^n ff-o(kms-') rTtT 



1 


0.1955 


0.6322 


1.47 


282 + 11 


2.38 ± 0.02 


2 


0.3317 


0.5235 


1.15 


349 ± 24 


3.37 ± 0.22 


3 


0.3223 


0.5812 


1.03 


326 + 16 


3.26 ±0.13 


4 


0.1642 


0.3240 


1.61 


325 + 12 


4.81 ± 0.02 


5 


0.2405 


0.4700 


1.32 


318 ± 17 


2.60 + 0.03 


6 


0.1260 


0.5349 


1.00 


229 ± 13 


1.82 + 0.05 


7 


0.2318 


0.7950 


1.15 


274 + 15 


1.77 + 0.01 


8 


0.0808 


0.7115 


0.85 


195 + 10 


1.23 + 0.01 


9 


0.2046 


0.4814 


1.39 


290 ± 16 


3.14 + 0.02 


10 


0.0629 


0.5352 


1.04 


206 ±5 


2.60 + 0.10 


11 


0.2076 


0.5241 


1.21 


295 ± 13 


2.14 + 0.02 


12 


0.2479 


0.7933 


1.81 


279 + 17 


2.02 + 0.02 


13 


0.2285 


0.4635 


1.25 


305 ± 19 


1.80 + 0.01 


14 


0.1553 


0.5170 


1.64 


271 ± 16 


4.20 + 0.04 


15 


0.0819 


0.5324 


1.57 


245 ±7 


4.47 + 0.01 


16 


0.485 


3.595 


1.34 


229 + 15 


0.82 + 0.12 


17 


0.938 


2.941 


1.24 


251 ± 19 


1.60 + 0.15 


18 


0.810 


3.399 


1.41 


224 ± 15 


1.06 + 0.08 


19 


0.497 


2.092 


0.36 


116 ± 10 


0.41 + 0.04 


20 


1.004 


3.263 


1.56 


328 ± 32 


0.31+0.06 



Refere nces - Treu et a l i2002] |. HQOU, ( l2004t . ( |2006] |: Koopmans et 
al. GmH . ( I2003ab , ( l2006l l. 



SLACS sample 



SLACS+LSD sample 




Fig. 7. Probability density function (/), marginalized over the 
values of the cosmological parameters, of the exponent (y) 
which characterizes the total density profile of the lens galaxies. 
The probability distribution with the approximate normal distri- 
bution for the SLACS (on the left) and for the SLACS+LSD (on 
the right) samples are obtained through MCMC methods. 



lens galaxies at redshifts zi ranging from 0.06 to 0.33 and back- 
ground sources with a redshift range of z, from 0.32 to 0.80. 
Detailed gravitational lensing models on the HST/ACS images 
were developed in order to measure the total mass inside the 
Einstein angle, i.e. 6e for an isothermal model, to a less than 
a few percent accuracy; from the same images also the values 
of the galaxy effective radii (Rg) were derived. Furthermore, 
(To and its uncertainty (average value of 5.0%) were obtained 
from the Sloan Digital Sky Survey Database. Afterwards, we 
studied a second data set, adding to the previous one five ad- 
ditional field elliptical galaxies from the LSD (Lenses Structure 
and Dynamics) Survey (Koopmans & Treu 120021 I2003al Treu 
& Koopmans |2002l |2U0^ |2TO] l (see Table©. These lens and 
source galaxies are at higher redshifts (zi ~ 0.5 - 1.0, Zs ~ 
2.0 - 3.6), thus it follows from Sect. 2 that the dependence of 
r on the cosmological parameters should be easier to test. On 
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Table 7. Mean value and standard deviation of the parameter y 
for the SLACS and SLACS+LSD samples. 



Cumulative distribution (ii„+fi^=l) 



SLACS 



SLACS+LSD 



1.99±0.04 



2.01±0.03 



the other hand, these galaxies have higher central velocity dis- 
persion errors, so that the mean uncertainty on o-q rises now to 
5.7%. 

At first, we checked on the two samples the plausibility of 
the hypothesis of homologous total density distribution, without 
any assumptions on the values of the cosmological parameters. 
In order to do so, we took into account the following relation 
(Koopmans |2005l l 



:; T = r(zi, z/, n,„, Qa) — - 



4-n 



Re 



2-y 



y(3 - 7) 



(4) 



which extends Eq. (|2|, considering a more general power law 
model for the total density profile (p oc 1/r^) of the lens galaxies 
(the isothermal case is retrieved setting -y equal to 2). The pa- 
rameter y was assumed to be the same for all the lenses and we 
used Markov chain Monte Carlo (MCMC) methods (with 5x10^ 
steps for each chain) to sample its probability density distribu- 
tion. The marginalized probability density function (/) with the 
approximate normal distribution, the mean value, and the stan- 
dard deviation of y for the two SLACS and SLACS-hLSD sam- 
ples are shown in Fig.|2]and Table|2] From these results, we can 
state that the total density profile of the lens ellipticals is indeed 
well approximated by an isothermal distribution (y equal to 2), 
independentely on the adopted cosmological model. 

Starting from here, we evaluated the observational angular 
diameter distance ratio of Eq. (|2]) and its uncertainty, assuming 
a reasonable 5% error on Oe- Hence, we used the x^ statistics 
reported in Eq. (O to gain information about the cosmological 
parameters. 

The results are plotted in Figs. [8]and|9] Here we show the 
intervals at the 95% confidence level for the dark-energy den- 
sity parameter in a flat cosmological model, and the 68%, 95%, 
and 99% confidence regions in the parameter space (Q„,, Qa) 
for general cosmological models. In a flat geometry, at 95% CL, 
Qa is found to be smaller than 0.80 from the SLACS sample, 
and greater than 0.16 from the SLACS+LSD sample. In addi- 
tion, without any assumptions on the cosmological parameters, 
the SLACS+LSD sample rules out, at greater than 99% CL, cos- 
mological models with dark-energy density parameter smaller 
than 0.4 (as shown in the bottom right panel of Fig. ^. From a 
comparison of Fig.[8]and Fig.|9l it turns out that the larger set of 
galaxies with more distant lenses shifts the minimum ;\f^ towards 
higher values of Qa and that a much bigger region of values for 
the cosmological parameters is excluded, at the same confidence 
level. 

The small number of lenses in the two samples is the primary 
reason that prevents us from reaching higher precision results. 
We note that the current concordance measurement of Qa sug- 
gests a value between 0.7 and 0.8 and this range is included in 
the 95% confidence intervals and in the 99% confidence regions 
for both the SLACS and the SLACS+LSD samples. 

6. Conclusions 

Lensing studies have been recognized to be very valuable in pro- 
viding a testing ground for theories of formation and evolution 
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Fig. 8. Results from the SLACS sample. Top left: Cumulative 
distribution (F) of the distance between the x^ in the true (Qa = 
0.7) and the best flat cosmological model from 2000 simula- 
tions with 15 lenses each. The dashed line gives the 95% con- 
fidence level. Top right: The x^ curve for the SLACS sample in 
a flat cosmology. The dashed line corresponds to the value of 
Ax^ of the previous panel and the thick bar on the abscissa axis 
shows the 95% confidence level interval for Qa. Bottom left: 
Cumulative distribution (F) of the distance between the x^ in 
the true (Q,„, Qa) = (0.3, 0.7) and the best cosmological model 
from 2000 simulations with 15 lenses each. The dashed-dotted, 
dashed, and dotted lines give, respectively, the 68%, 95%, and 
99% confidence levels. Bottom right: The x^ contour plot for 
the SLACS sample; the levels correspond to the values of lS.x^ of 
the previous panel. The dashed-dotted, dashed, and dotted lines 
show, respectively, the 68%, 95%, and 99% confidence regions 
for Qm and Qa. 



of early-type galaxies (Rusin & Kochanek 120031 Koopmans et 
al. I2OO6I 1. In these studies the structure of the lens galaxies can 
be investigated, once a specific cosmological model is adopted. 
From a complementary perspective, a measurement of the cos- 
mological parameters can in principle be accomplished by as- 
suming some empirically justified general properties for a sam- 
ple of lenses. In this paper, we have shown that indeed, starting 
from the 1 jr^ paradigm for the total density profile of elliptical 
galaxies, the cosmological parameters Q„, and Qa can be esti- 
mated. The adopted paradigm has been validated by stellar dy- 
namics (e.g., Gerhard et al. 1200 11 1, strong gravitational lensing 
(e.g., Rusin et al.'2003), and dynamical studies based on several 
other tracers: globular clusters and planetary nebulae (e.g., Peng 
et al. |2004] |. X-ray halos (e.g., Humphrey et al. |2006] |. and HI 
disks and rings (e.g., Franx et al. .l994 ). 

In detail, the main results of this paper can be summarized 
as follows: 

- We have developed a new method to investigate the geom- 
etry of the Universe by combining measurements of strong 
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(isothermal) model for the total density distribution in ellip- 
ticals. For a sample of eight bright, round, and nearby early- 
type galaxies, we have shown that the stellar central velocity 
dispersion is a good tracer of the velocity dispersion charac- 
terizing the one-component model. 

- We have applied the proposed method to the SLACS Survey 
data set. This sample is composed of fifteen massive early- 
type galaxies at intermediate redshift acting as lenses. In a 
flat cosmology this method leads to a value of Oa lower than 
0.80 at 95% CL. By including five more distant lenses from 
the LSD Survey, at the same confidence level and in the same 
flatness hypothesis, Qa has been measured to be greater than 
0.16. This latter sample ruled out, at greater than 99% CL, 
general cosmological models with values of Qa smaller than 
0.4. 

- Finally, we conclude that future surveys expected to iden- 
tify a large number of lenses (one or two hundred) will al- 
low this method to measure the values of the cosmological 
parameters with a precision comparable to that of other stan- 
dard techniques. This will be a new and important test for the 
KCDM concordance model. 
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Fig. 9. Results from the SLACS+LSD sample. Top left: 
Cumulative distribution (F) of the distance between the x^ in 
the true (Qa - 0.7) and the best flat cosmological model 
from 2000 simulations with 20 lenses each. The dashed line 
gives the 95% confidence level. Top right: T\\t y^ curve for the 
SLACS+LSD sample in a flat cosmology. The dashed line cor- 
responds to the value of Ax^ of the previous panel and the thick 
bar shows the 95% confidence level interval for Qa- Bottom left: 
Cumulative distribution (F) of the distance between the x^ in 
the true (Q,„,Qa) - (0.3,0.7) and the best cosmological model 
from 2000 simulations with 20 lenses each. The dashed-dotted, 
dashed, and dotted lines give, respectively, the 68%, 95%, and 
99% confidence levels. Bottom right: The x^ contour plot for 
the SLACS+LSD sample; the levels correspond to the values of 
A;^f^ of the previous panel. The dashed-dotted, dashed, and dot- 
ted lines show, respectively, the 68%, 95%, and 99% confidence 
regions for Q.^ and Qa- 



gravitational lensing and stellar dynamics in a sample of el- 
liptical galaxies. The basic idea is to compare lensing and 
dynamical mass estimates in order to find an observable 
cosmology-dependent relation. For this purpose, it is nec- 
essary to model the total (luminous and dark) density profile 
inside the Einstein radius of the ellipticals under investiga- 
tion. Here, we have assumed a one-component 1 jp- density 
profile. 

We have performed a feasibility study of our technique by 
simulating measurements of the cosmological parameters. 
We have studied the precision of the method, which depends 
on the level of uncertainty with which the Einstein angle and 
the central stellar velocity dispersion are known. Hence, we 
have demonstrated that the current uncertainty with which 
these quantities can be measured is sufficiently small to ob- 
tain precise information about cosmology, once a statistically 
significant sample of lenses is available. In particular, the 
method has been shown to be best suited to measure Qa. 
Then, we have studied the relation between the stellar kine- 
matics and the velocity parameter of the one-component 1 /r^ 
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